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Increasing improvements in the independent determinations of the Hubble 
constant and the age of the universe now seem to indicate that we need 
a small non-vanishing cosmological constant to make the two independent 
observations consistent with each other. The cosmological constant can be 
physically interpreted as due to the vacuum energy of quantized fields. To 
make the cosmological observations consistent with each other we would need 
a vacuum energy density, p„ ~ {IQ^^eVY today ( in the cosmological units 
h = c = k = l). It is argued in this article that such a vacuum energy density 
is natural in the context of phase transitions linked to massive neutrinos. In 
fact, the neutrino masses required to provide the right vacuum energy scale 
to remove the age Vs Hubble constant discrepancy are consistent with those 
required to solve the solar neutrino problem by the MSW mechanism. 



1 Introduction 



Increasing accuracy in astronomical observations is leading us to an increasing precision 
in the determination of cosmological parameters. This in turn is leading us to critically 
re-examine our cosmological models. In particular, the precise determination of the 
Hubble constant and the independent determination of the age of the universe is forcing 
us to critically re-examine the simplest and most appealing cosmological model - a flat 
universe with a zero cosmological constant [|I], 0. 

The Hubble constant enters in the relationship between the recession velocity of an 
object and it's distance from us. The recession velocity of an object can be determined 
by using the Doppler effect and is relatively easy to determine. It is the calibration 
of the extragalactic distance ladder which is the difficult part of measuring the Hubble 
constant and in which the precision has been increasing significantly. It is the use of 
the Cepheid variables as standard candles that has allowed the improved determination 
of the extragalactic distance scale. Cepheids are variable stars whose pulsation period 
are very stongly correlated with their luminosities. These stars are well understood 
theoretically and the period - luminosity relationship is well-documented empirically. 
By observationally determining the pulsation period of a Cepheid variable and using 
the period luminosity relationship one can immediately determine the luminosity of the 
object. Then by using the apparent brightness of the object one can accurately determine 
the distance to the Cepheid. Pierce et al|l|] have used this technique with ground-based 
observations to determine the extragalactic distance scale. An excellent discussion on 
this subject is contained in the article by Pierce et al^ and references therein. 

Pierce et al|l]] have determined the Hubble constant to be, = 87 ±7 km s~~^ Mpc~^. 
They further point out in their article that this value of the Hubble constant is in fact in 
conflict with the independent determination of the age of the universe]^] using Galactic 
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globular clusters if we use the standard cosmological model with a zero cosmological 
constant. The estimate of the age derived from an analysis of the galactic globular 
clusters is 16.5 ± 2 Gyr. 

An accurate determination of the Hubble constant is also one of the important goals 
of the Hubble Space Telescope (HST). In fact, Freedman et al have used the HST 
to calibrate the extragalactic distance scale and hence determine the Hubble constant. 
They obtainn the value of the Hubble constant to be 80 ± 17km s^^ Mpc^^. They also 
point out that their determination of the Hubble constant is inconsistent with the age 
of the globular clusters within the framework of standard f2 = 1 cosmology with no 
cosmological constant. 

2 Resolution of the age Vs. Hubble constant prob- 
lem through the introduction of a small vacuum 
energy. 

One of the ways to avoid the apparent conflict between the observed age of the uni- 
verse and the observed Hubble constant is to introduce a small cosmological constant 
in the Einstein equations that govern the evolution of the universe. This idea has been 
extensively studied by a number of people including Tayler|^ and Klapdor and Grotz[||. 

Let's quickly summarize how a cosmological constant of the right magnitude can solve 
the apprent conflict between the age and the Hubble constant. This can be seen from the 
following analysis^]. To a good approximation our universe is spatially homogenous and 
isotropic on large scales. It is therefore appropriate to describe space-time by Robertson- 
Walker metric which can be written in the form, ( units c = 1 ) 



ds' = dt'-R\t) 



r^d9^ + sin^ 



1) 



— kr'' 

where, (t, r, 9, 0) are the comoving coordinates describing a space time point and 
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R{t) is the cosmic scale factor. Also, k = +1, —1 or depending on when the universe is 
closed, open or fiat. 

The present expansion age of a matter dominated universe can be evaluated in a 
Robertson- Walker universe. General fine-tuning arguments as well as the inflationary 
picture gives us a preference for a fiat universe, = 1. In this case, to = ^H~^. For, 
fig ~ 1, one can expand the above expressions in a Taylor expansion. 



to = Ih;' 



i-^(^^o-i) 
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We can also determine the present age of the universe containing both matter and 
vacuum energy such that Q^ac + ^matter = 1, 



to = Ih;'^;;^/' in 



1 + Q^/"^ 



(3) 



^1 fi^QIf.)"'^^^ 

This will give us much longer lifetimes as can be seen most dramatically by examining 
the limit Qvac 1 in which case to — > oo. Indeed having an Qvac ~ 0.8 is one of the 
possible solutions of the age Vs. Hubble constant discrepancy as can be seen from the 
following discussion. 

First, let us quicldy recall the observational numbers on the Hubble constant and the 
ages of globular clusters. Here it is worthwile pointing out that there are actually 2 sets 
of numbers which though consistent with each other have slightly different central values 
and error estimates. Pierce et al[|| quote the result of their analysis as yielding a Hubble 
constant. Ho of 87 ± 7km Mpc~^. They then draw attention to the fact that this 
is in conflict with the age estimate of the globular clusters inferred by VandenBergh'sP| 
analysis which gives an age of 16.5 ± 2 Gyr. 

Freedman et alp| have a greater amount of data that they have analyzed very thor- 
oughly. Thus they found over 20 Cepheids in a Virgo cluster galaxy as opposed to 3 
Cepheids found by Pierce et al. They obtain a slightly different central value of Ho with 



substantially larger and perhaps more realistic error bars than Pierce et al. They thus 
quote a Hubble constant value, Ho = 80 ± llkm Mpc~^. Also based on a wider 
spectrum of data for ages they quote a central value of the age of the universe as 14 Gyr 
with 1 (7 error bars of ±2 Gyr. Preedman et al, however, too point out that even with 
their more generous error bars there is still a discrepancy between the Hubble constant 
and the age of the universe if we restrict ourselves to a standard fl — 1 cosmology with 
a zero cosmological constant. 

Let us now see how a non-zero cosmological constant can solve this discrepancy and 
what values of the cosmological constant are typically required to minimally solve the 
discrepancy between the Hubble constant and the age of the universe. First, let us 
consider the lowest possible value of the Hubble constant as quoted by Freedman et al 
which is 63 km Mpc~^. Let us consider this together with the age of the universe, 
to — 14.5 Gyr which though slightly higher than the central value quoted by Freedman 
et al is at the lower hmit of the age quoted by Pierce et al. These two values of Ho 
and to imply a value of Ho to = 0.93. This corresponds to the value of flyac = 0.66. 
Thus, Q^ttc = 0.66 would remove the contradiction between the Hubble constant as 
determined by Freedman et al and the age of globular clusters. However, let us consider 
a few more values of flvac, Ho and to to see where future more precise observations 
and analysis might lead us. Consider Q^ac — 0.7 which imphes Ho to — 0.964. For 
Ho — 63km Mpc~^ this would imply the age to = 15 Gyr and for to = 14.5 Gyr this 
would imply a Ho = 65km Mpc~^. Similarly, flyac — 0.8 would imply Ho to = 1.076 
which would give an Ho — 73km s"^ Mpc~^ for to — 14.5 Gyr and for a, to — 16.5 Gyr 
would imply Ho = 64A;m Mpc'^. Finally, if we consider the lower bound implied by 
the numbers quoted by Pierce et al, viz. Ho = 80km Mpc~^ and to = 14.5 Gyr, we get 
^vac ~ 0.85. Clearly however one would not like to push the values of Q^^c much higher 
than this number for a number of reasons. First, we would like fl^ac + ^matter — 1 and 
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higher values of Qvac will start to conflict with the lower bound on matter density from 
galaxies and clusters. Furthermore, one can start to place an independent constraint 
on the cosmological constant fromm gravitational lens statistics. Thus based on the 
HST Snapshot survey of quasars Maoz and Rix|S^ have inferred the bound Q^ac < 0.7. 



Fukugita, Hogan and Peebles 1 36] have stated that this constraint may be relaxed slightly 
if one considers the fact that at large distances typical of lensing events galaxies are 
observed to have larger star formation rates than their nearby counterparts and may also 
contain larger quantities of dust. This may lead to some lensing events getting obscured 
and hence the upper limit for Q^ac niay rise if this is taken into account. It is unclear 
how much the upper limit will be pushed upwards but the fact remains that for both the 
above reasons one would not like an extremely large value of ^l^ac even though a value of 
Q^ac ~ 0.8 perhaps best meets all the observational constraints outlined above. Improved 
more accurate results and more exhaustive analysis on all these fronts will clearly shed 
further light on all these connected issues. Let us now turn our attention to a discussion 
of the possible physical origin of such an flvac- 
The value f2 

yac ~ 0.8 corresponds to an energy density of the vacuum energy density 
of py ~ This energy density is of course much lower than most familiar 

energy scales in particle physics and the question naturally arises as to the physical origin 
of this energy scale. The smallness of this energy scale has been frequently referred to as 
the cosmological constant problem. 

However, we'll argue in this paper that in fact a cosmological constant of the right 
magnitude required to make the cosmological observations consistent with each other 
may follow from the dynamical evolution of our universe. The basic physical picture 
which will allow us to arrive to this conclusion is that the cosmological constant might 
be interpreted as the vacuum energy of the quantized fields. 

This point has been made by many people and is discussed at length by Birell and 



Davies]^. Further, this vacuum energy is not a static quantity but a function of time. 
This idea too has been extensively explored by a number of people including Peebles and 
RatralBl, Freese, Adams, Frieman and MottolalHl, Renter and Wetterich||lU||. In fact, we 



know that there were a number of phase transitions in the evolution of the universe. 

Thus, the history of the universe may be summarized as periods of dramatic change 
charecterized usually by phase transitions with relatively quiet periods of relaxation 
between the phase transitions [Q. Indeed, since the vacuum energy density changes 
as the {characteristic energy scale)^ at a phase transition, in the absence of fine- 
tuning one expects that the vacuum energy density at the end of a phase transition 
~ {characteristic energy scale)^ 

This idea has been spelt in detail in the paper by Wilczek[0] and also by Renter 
and Weterich||10[. Thus, p^, ~ (10^^ GeV)"^ at the Grand Unified Symmetry breaking, 
Pv ~ (10^ GeVY at the elctroweak symmetry breaking and ~ (10^^ GeVY at the 
chiral symmetry breaking in QCD. 

Furthermore, at the conclusion of a phase transition the vacuum energy starts decay- 
ing more slowly to the energy scale characterized by next phase transition. This point of 
view is implicit in the papers by Wilczek|^l| and has been explicitly stated by Renter and 



Wetterich. In fact, the physical mechanism for the decay of vacuum energy is coupling 
to lighter fields. This mechanism is briefiy discussed by Freese, Adams, Frieman and 
Mottolap, who do an extensive analysis of cosmology with a decaying vacuum energy. 

Thus, the question of the magnitude of the cosmological constant really becomes a 
question about energy scales. Almost every paper on the subject of the cosmological 
constant has had to struggle with the characteristic energy scale of ~ lQ~^eV in the form 
p^ ~ (O.OOSey)^ in various guises such as (IQ-^^Gel/^) or A/M^^ ~ lO^^^o 

The fact that O.OOSeV is so much less than any characteristic energy which familiar 
to most of us has caused a great deal of consternation. 



However, the energy scale 0. 003eV is certainly not a complete stranger to us. The 
most natural low-energy scale that particle physics gives us is the light neutrino masses 
that follow from the see-saw model of neutrino masses. In fact, the neutrino masses 
required to to solve the solar neutrino problem by MSW mechanism[^] imply neutrino 

masses ~ 10~^eV 

This, of course, is a powerful hint but is not yet a solution to the cosmological constant 
problem. 

In fact, the finite temperature behaviour of the see-saw model of neutrino masses has 
been studied in detail by Holman and Singh |26|]. The original motivation for studying 



this model was to provide a concrete particle physics model fo the Late Time Phae 
Transitions model for structure formation. Our analysis showed that in fact this model 
does exhibit a phase transition with a critical temperature ~ {few) m^. 

3 Late Time Phase Transitions and the time evolu- 
tion of the Hubble parameter 

In this section we will discuss the cosmological motivations and particle physics models 
for Late Time Phase Transitions. Once we have a specific model we will study the time 
evolution of fields and the scale factor in this model. In particular, we will be interested 
in studying the time evolution of the Hubble parameter and will see that in this model the 
Hubble " constant ", in fact has an acceptable value at the present age of the universe. 

Phase transitions that occur after the decoupling of matter and radiation have been 
discussed in the literature as Late Time Phase Transitions (LTPT's). The original mo- 



tivation for considering LTPTsM] [|T9[ pOl p^] [E^ [E6|was the need to reconcile the 



extreme isotropy of the Cosmic Microwave Background Radiation (CMBR)|[T5| with the 
existence of large scale structure [^] and also the existence of quasars at high redshifts|r?| . 
Discussions of realistic particle physics models capable of generating LTPT's have 



been carried out by several authors p3| [p2| . It has been pointed out that the most 
natural class of models in which to realise the idea of LTPT's are models of neutrino 
masses with Pseudo Nambu Goldstone Bosons (PNGB's). The reason for this is that the 
mass scales associated with such models can be related to the neutrino masses, while any 
tuning that needs to be done is protected from radiative corrections by the symmetry 
that gave rise to the Nambu- Goldstone modes [^ . 

Holman and Singh[^ studied the finite temperature behaviour of the see-saw model 
of neutrino masses and found phase transitions in this model which result in the formation 
of topological defects. In fact, the critical temperature in this model is naturally linked 
to the neutrino masses. 

The original motivation for studying the finite temperature behaviour of the see-saw 
model of neutrino masses came from a desire to find realistic particle physics models for 
Late Time Phase Transitions. It now appears that this may also provide a physically 
appealing and observationally desirable magnitude for the cosmological constant. 

In particle physics one of the standard ways of generating neutrino masses has been 
the see-saw mechanism [^. These models involve leptons and Higgs fields interacting 
by a Yukawa type interaction. We computed the finite temperature effective potential of 
the Higgs fields in this model. An examination of the manifold of degenerate vacua at 
different temperatures allowed us to describe the phase transition and the nature of the 
topological defects formed. 

To investigate in detail the finite temperature behaviour of the see-saw model we 
selected a very specific and extremely simplified version of the general see-saw model. 
However, we expect some of the qualitative features displayed by our specific simplified 
model to be at least as rich as those present in more complicated models. 
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3.1 A particle physics model for LTPT 

We chose to study the 2-family neutrino modeL Because of the mass hierarchy and small 
neutrino mixings [jl2l we hope to capture some of the essential physics of the Ve-^n system 
in this way. The 2-family see-saw model we consider requires 2 right handed neutrinos 
which transform as the fundamental of a global SUr{2) symmetry. This symmetry 
is implemented in the right handed Majorana mass term by the introduction of a Higgs 
field o"jj, transforming as a symmetric rank 2 tensor under 5'f//j(2) (both and aij are 
singlets under the standard model gauge group). The spontaneous breaking of SUr{2) via 
the vacuum expectation value (VEV) of a gives rise to the large right handed Majorana 
masses required for the see-saw mechanism to work. Also, the spontaneous breaking 
of SUr{2) to U{1) gives rise to 2 Nambu Goldstone Bosons. The SUr{2) symmetry is 
explicitly broken in the Dirac sector of the neutrino mass matrix, since the standard 
lepton doublets II and the Higgs doublet $ are singlets under SUr{2). It is this explicit 
breaking that gives rise to the potential for the Nambu Goldstone modes via radiative 
corrections due to fermion loops. Thus, these modes become Pseudo Nambu Goldstone 
Bosons (PNGB's). 

The relevant Yukawa couplings in the leptonic sector are: 

- C,^^ = ya^'iN'j,^ + yNlNi V,, + h.c. (4) 

where a,2,j = 1,2. The SUii{2) symmetry is implemented as follows: 

N'r - f/jiv^ 

- u^Mu^y, (5) 

where t/j is an SUfi{2) matrix. The first (Dirac) term above transforms as an SUr{2) 
doublet, thus breaking the symmetry explicitly. 
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We now choose the VEV of o to take the form|]2^: ioij) = f6ij, thus breaking SUr{2) 
spontaneously down to the U{1) generated by T2 (where Xj are the Pauh matrices). We 
take / to be much larger than the doublet expectation value v. 

We can parametrize cxjj so as to exhibit the Nambu-Goldstone modes as follows: 

a{x) = U{x){a)U^{x) 

= fU{x)U^{x) (6) 

with U{x) = exp(i(^iri + ^3X3)//) (note that U is symmetric). 

After the Higgs doublet acquires its VEV, we have the following mass terms for the 
neutrino fields: 

- £^ass = rUaiulNi, + MN rUU^ N% + h.c. (7) 

where are the standard neutrinos, ruai = Vai f^/v^, M = yf/\/2. 

Diagonalizing the neutrino mass matrix in the standard see-saw approximation (|m(jj| << 
M) and performing a chiral rotation to eliminate the 75 terms, we find that the de- 
pendent light neutrino masses are given by 

m\ = ^[cos22||^|| (m2,+my2 + sin22||e|| (f3(m2,-m22) + 2mnmi2fi)'] 
ml = ^[cos2 2||e|| (ml^ + mj^f + sin^ 2m\ (^(m^, - m^s) + 2m2im22fi)'] 

(8) 

where ||^|| = ^J^f + ^3// and = ii/\Jii + ^3. We neglect the effects of the heavier 
neutrinos since they will be suppressed by powers of m^j /M in the loops that will generate 
the effective potential for the ^j's. 

For simplicity, to begin with we shall restrict ourselves to the case where the Dirac 
mass matrix mai is proportional to the identity: m^j = m 5ai- We will consider a more 
general case later. Using standard results on the computation of the effective potential 

10 



due to fermion loops[Q (we treat the ^j's as classical background fields, i. e. we do not 
allow them to propagate in loops), we can calculate the one- loop effective potential for 
the ^j's. The renormalized potential can be expressed as follows: 



^re„orm(ei, ^s) = K) + m^M' + KiM')' - — (^^)^(ln ^ - -) (9) 



where /i is the subtraction point. Note that we can absorb the effects of A,, by redefining 
fi. We will suppose that this has been done in what follows. Further, the quantity A4^, 
is given by 



-M' = ^(cos^2||e||+gsin2 2||e||). (10) 



The finite temperature correction due to the two light neutrinos is given by, 



T^4 »oo 

AV^r(ei,e3) = -4— / c^xx^log 



1 + exp 



2 



1/2 



(11) 



2^2 

The above expression can be evaluated numerically for any given choice of parameters, 
however for some purposes it is useful to expand the above expression to get analytic 
expressions. Performing the high temperature expansion of the complete potential and 
discarding terms of order (A^^)^/T^ or higher we get. 



VUCi, ^3) = V{M') = [Vo - — — ) + {ml + Ty6)M' + ^-^{n - log -),(12) 

90 Stt'^ 



where n = 27 — 1 — 21og7r ~ —2.1303, rrir is a parameter in the model and /i is the 
renormalisation scale. A4 is naturally of the neutrino mass scale in this model. 

A study of the manifold of degenerate vacua of the effective potential at different 
temperatures revealed phase transitions in this model accompanied by the formation of 
topological defects at a temperature of a few times the relevant neutrino mass . Typically 
at higher temperatures the manifold of degenerate vacua consisted of a set of disconnected 
points whereas at lower temperatures the manifold was a set of connected circles. Thus, 



11 



domain walls would form at higher temperatures which would evolve into cosmic strings 
at lower temperatures. 

Since the critical temperature of the phase transition, ~ (few) , let us quickly 
summarise the observational evidence for a small m^. 

At neutrino detectors around the world, fewer electron neutrinos are received from 
the sun than predicted by the Standard Solar Model. An explanation of the deficiency is 



offered by the MSW mechanism ||12|| which allows the Ve produced in solar nuclear reac- 
tions to change into v^. This phenomenon of neutrino mixing requires massive neutrinos 



with the masses for the different generations different from each other |12 



The model we considered earlier was an extremely simple one. Although it had 2 
families of light neutrinos, there was only one single light neutrino mass. As such this 
model was not compatible with the MSW effect. However it is fairly straightforward to 
modify our original model to make it compatible with the MSW effect as is shown in 
what follows. 

To ensure that it is not possible to choose the weak interaction eigenstates to coincide 
with the mass eigenstates we must require the 2 neutrino mass scales to be different. We 
can ensure neutrino mixing in our model by demanding that irtai be such that mn ^ 17122 
and mi2 = = 777,21. In this case, the effective potential KotlG^^s) = l/2(F(A^f) + 
V{Ml)) with V{M'^) having the same functional form as V{M'^){i = 1,2) and M'^ 
given by the following expression: 

M^ = ^(cos^2||e||+gsin^2||e||) (13) 



Further, if mn << 17122 then Vtot(^i,^3) = V{-M.l)/2, which is exactly half the finite 
temperature effective potential we discussed earlier except the neutrino mass scale is the 
heavier neutrino mass scale. Hence, the discussion on phase transitions and formation 
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of topological defects we carried out earlier goes through exactly except that the critical 
temperature is determined by the mass scale of the heavier of the 2 neutrinos. 

In the complete picture of neutrino masses the neutrinos might have a mass 
hierarchy analogous to those of other fermions. Further, we expect that the mixing 
between the first and third generation might be particularly small . In this scheme, it is a 
good first approximation to consider 2-family mixing. We are here particularly interested 
in the Ve-^fi mixing. This is also the mixing to which the solar neutrino experiments are 
most sensitive. A complete exploration of MSW solutions to the solar neutrino problem 



has recently been reported by Shi, Schramm and Bahcall[^. We shall restrict ourselves 
to the 2-family mixing. The data seems to imply a central value for the mass of the muon 
neutrino to be a few meVpOj. 

We now turn to a quick discussion of the distortions of the Cosmic Microwave Back- 
ground Radiation (CMBR) this model produces. The most significant microwave distor- 
tion comes from collapsing domain wall bubbles. This has been discussed and calculated 
by Turner, Watkins and Widrow (TWW) [|T^] . As pointed out by TWW this anisotropy is 
most significant on ~ 1° angular scales. The temperature shift due to a photon traversing 
a collapsing domain wall bubble is 
AT 

— = 2.64 X 10~^/i-^/?AcT/(10Me1/=^) (14) 

where h. A, (3 are dimensionless numerical constants of order unity and a is the surface 
tension of the domain wall. The present measurements of the CMBR anisotropy then 
implyg that a < O.SMeK^. 

An estimate of a in terms of the quantities mi, and / introduced in our model can be 
obtained|T^. (To make contact with the work of L. Widrow cited above please note that 
his \m* = m{v^) and m = / in our notation. ) Thus, the constraint on a then implies 
that / < l{]^^GeV . Our model is clearly an effective theory with / being some higher 
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symmetry breaking scale on which it is tough to get an experimental handle. However, 
the constraint derived above is in fact natural in the context of the see-saw model of 
neutrino masses embedded in Grand Unified Theories as discussed by Mohapatra and 
Parida (MP)[gg] and also by Deshpande, Keith and Pal(DKP)|3§. 

3.2 Time evolution in the LTPT model 

Now that we know the potential in which the fields and ^3 evolve, we can write down 
the coupled set of evolution equations which describe the time evolution of the fields and 
the scale factor of the universe. Once again we'll follow the general techniques described 
in Kolb and Turner [Q. The time evolution of the scale factor is given by equations like 
(2), (3) and (4). It is perhaps worth noting that the expression for the presssure and 
energy density of the fields is given by: 

p« = ^^^ + nei,e3) (15) 

Pi = ^-^^-V{^u^s) (16) 
The time evolution of the fields .^1 and .^3 in an FRW universe is given by. 

These coupled equations describing the time evolution can be solved numerically. Here 
we are interested in the time evolution at very recent epochs. Clearly, for extremely recent 
epochs the high temperature expansion is inappropriate. However, the zero temperature 
potential which we have computed and described is a good approximation for stuudying 
the time evolution at recent epochs. Thus we will use the zero temperature potential to 
do the time evolution in what follows. In fact as it turns out the time evolution of the 
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scale factor is fairly insensitive to the initial conditions on the fields but is determined 
primarily by the order of magnitude of the energy density in the fields. In fact we 
evolved the system with a variety of initial conditions on the fields and observed an 
almost identical time evolution for the scale factor. 

We have investigated the behaviour of the system for a variety of choice of parameters 
entering into the potential with the typical order of magnitude for the parameters ~ 
nil, and Vg ~ m^. Since it is the fact that the vacuum energy density in our model is 
~ mf, that plays a crucial role in the time evolution in our model, let us re-emphasise 
why this is natural in the context of our model. 

Recall, that the original — >Cyuk contained 2 distinct couplings, i/ai and y. In particular, 
if we set ijai = in the Lagrangian then the symmetry of the Lagrangian was enhanced, 
the light neutrino masses would have been identically zero and the potential for the 
fields would have also been identically zero. If yai 7^ then the SUr{2) symmetry is 
explicitly broken, the light neutrino fields pick up their masses as outlined earlier and also 
the ^i's develop a non-trivial potential. Furthermore, since setting yai — enhances the 
symmetry of the Lagrangian, there is a symmetry which protects the small parameters 
in this model. In the model we are studying the coupling of the light neutrino fields to 
the ^ fields is identically zero if yai — and arises at the second order in perturbation 
theory in the see-saw diagonalization if yai 7^ 0. It is the coupling of the neutrino fields 
to the ^ fields that is responsible for a non-zero effective potential for the $, fields, hence 
there is a prefactor to the entire effective potential (including the vacuum energy part) 
which is proportional to the appropriate power of yai- The observable quantity that 
yai corresponds to is rriiy as given in section 3.1 .Since the effective potential has four 
mass dimensions, the dimensional prefactor multiplying the potential ends up being m^. 
Thus it is natural that the contribution to the vacuum energy density due to the fields 
appearing in the see-saw model of neutrino masses presented here ends up being ~ mf,. 
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One, of course, has to worry about the contribution of other heavier fields to the vacuum 
energy density in the cosmo logical context. The underlying picture being used is that 
discussed by Wilczek|TT]], Renter and Wetterich[^] and Freese, Adams, Frieman and 



Mottola[^. They argue that the cosmological constant will decay during the evolution 
of the universe as the vacuum energy of heavier fields dissipates due to their coupling 
to lighter fields. One would then expect the vacuum energy density at late times to be 
dominated by the contribution of the lightest most weakly coupled fields such as those 
appearing in the see-saw model of neutrino masses discussed in this paper. An in-depth 
analysis of the details of this mechanism will be the subject of a later work. 

The time evolution of the system can be summarized as follows. The fields evolved 
to the minimum of their potential on a time scale which is short compared to the typical 
Hubble time scale in the problem. The evolution of the scale factor follows the normal 
matter dominated behaviour for a while until the vacuum energy starts playing an im- 
portant role. After this time the vacuum energy starts driving the time evolution of 
the scale factor. Thus it is the value of the vacuum energy density that determines the 
asymptotic time evolution of the system. For our model with the choice of parameters 
stated and rationalised above we have the final vacuum energy density ~ m^. 

The observationally important plot is the plot of the Hubble parameter as function 
of time. This is displayed in figure 1. As you can see the Hubble parameter assumes a 
constant value after the vacuum energy starts playing the dominant role in the evolution 
of the scale factor. 

In fact, in retrospect one can understand the time evolution of the coupled differential 
equations simply by noting the order of magnitude of the quantities involved in the 
evolution equations. 

Let us introduce the following dimensionless physical quantities, t = H t , rji = ^ and 
V = We'll use the following physical quantities to make the dimensionless quantities 
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of order 1, H ^ 75 h km sec'^ Mpc-\ f = lO^^ /12 GeV and = 2.5 ma x IQ-^eV. 

Here are some of the quantities of interest and their magnitudes. First, the expressions 
for the pressure density and energy density of the fields is given below: 



Pr)^rnl< 6.5 x 10 



m' 



p, = mt{ 6.5 X 10-''^ 

171 J 



dr I \ dr I 



+ ^(771,773) 



(19) 



(20) 



Since we have scaled quantities so that , h, r]i, t and V are all of order 1 it follows 
that, 



Pe~-y(ei,6) 



(21) 
(22) 



Thus in fact, for all practical purposes, we have 



(23) 



which is the equation of state for vacuum energy and hence demonstrates that this 
solution is very close in spirit to the cosmological constant solution for the age vs Hubble 
constant problem. What we have achieved is to provide a physical basis for the correct 
order of magnitude for this effective cosmological constant. This can be seen clearly by 
making the evolution equation for the scale factor dimensionless too. 



+ 



k 



1 dR{T) 
R{t) dr ' R^{t)H^ 



2 f2 
12 



'dr]i 
dr 



+ 



'drjs 
dr 



+ 0.85^V(77i, 773) 
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Note that it is the fact that the mass of the neutrinos is the correct order of magnitude 
which allowed the age of the universe to become compatible with the Hubble constant 
observed today. 

The picture presented in this article is of a vacuum energy that changes as a function 
of time due to the coupling of the fields responsible for the vacuum energy to the other 
fields. One may worry therefore that the vacuum energy may disappear because of the 
coupling of the ^ fields to other fields. However because ^ can only decay into the lighter 
neutrinos it is coupled to, the time scale on which this vacuum energy will dacay is much 
larger than the present age of the universe. This can be seen by calculating the decay 
width of the ^ , T^. 

This decay width of the ^ particles arises because of the coupling of the ^ to the 
lighter fermions with coupling y is given byQ] 

r, = ^ (24) 

where, 

ml^^ (25) 
Therefore, 

mj^'f (26) 

Thus the timescale on which the energy in the ^ fields is converted into the energy of z/'s 
is given by 

rr - (27) 

^ y mi/ nil, 

which is much greater than the present age of the universe. Thus this vacuum energy is 
clearly not a short lived thing. As the expression above displays this is a consequence 
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both of the fact that we have hght particles involved and that they are extremely weakly 
coupled to other particles. 

In conclusion, the MSW solution to the solar neutrino problem seems to imply a muon 
neutrino mass of a few meV. This in turn would lead to a phase transition in the PNGB 
fields associated with massive neutrinos with a critical temperature of several meV. This 
would then give us a vacuum energy density ~ (10~^ eV)^, which would help resolve the 
conflict between the independent determinations of the Hubble constant and the age of 
the universe. This phase transition also happens at the correct epoch in the evolution 
of the universe to provide a possible explanation of the peak in quasar space density at 
redshifts of 2 to 3[|T . 
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FIGURE CAPTIONS 

Figure 1 : Time evolution of the Hubble Parameter in the LTPT model : 
Time. 
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